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I.  INTRODUCTION 


Many  calculations  of  transient  solid  propellant  regression  rates 
assume  a  quasi-steady  dependence  on  pressure  (and  perhaps  initial  temp¬ 
erature).  Constant  volume  bombs  and  constant  pressure  strand  burner 
measurements  are  the  usual  source  of  the  supporting  data.  Experiments, 
however,  have  demonstrated  that  in  a  rapidly  changing  pressure  field, 
the  regression  rate  can  deviate  substantially  from  the  quasi-steady  rate 
for  increasingly  Qr  decreasing3>4  or  oscillating  pressure.  Because  the 
immediate  interest  is  in  increasing  pressures  found  in  gun  chambers,  the 
body  of  theory  and  experiment  for  oscillating  pressures  will  be  largely 
ignored. 

As  observed  in  Kuo*s^  recent  review,  models  which  rely  only  on 
rate  of  change  of  pressure  are  valid  only  when  the  excursions  are  small 
enough  for  linear  analysis.  For  larger  excursions  from  steady  state, 
a  thermal  theory  model  has  been  frequently  used  wherein  only  the  gas 
phase  is  assumed  quasi-steady.  For  extreme  excursions  even  this  approxi¬ 
mation  must  be  abandoned  although  estimates  have  been  seriously  proposed 
of  the  application  limits  of  the  quasi-steady  gas. 

Nelson^  has  shown  that  for  a  given  gun  pressurization  history, 
calculated  transient  rates  can  vary  widely  with  the  model  used.  Kooker 


1 

J.  Brulard>  P.  Kuentzmann^  R.  Kling >  n Reponse  drun  Propergol  Solide 
a  un  Echelon  de  Pression”y  La  Recherche  Aerospatiale E,  279-287  (1975). 

2 

S.  L.  Turkj  R.  A.  Battista >  K.  K.  Kuo>  L.  H.  Caveny and  M.  Summerfield, 
” Dynamic  Responses  of  Solid  Rockets  during  Rapid  Pressure  CHange”^  J. 
Spacecraft  &  Rockets >  10j  137-142  (1973). 

3 

C.  F.  Yin  and  C.  E.  Hermance ^  rfContinous  Measurement  Transient  Burning 
Rates  of  a  Composite  Propellant  Undergoing  Depressurization” ,  AIAA 
Paper  71-173  (1971). 

4 

C.  E.  Woolridge  and  G.  A.  Marxmany  ”A  Comparison  Between  Theoretical  and 
Experimental  Extinction  Behavior  of  Composite  Solid  Propellants” 3  AIAA 
Paper  70-666  (1970). 

5 K .  K.  Kuo  and  G.  R.  Coates >  ” Review  of  Dynamic  Burning  of  Solid  Propel¬ 
lants  in  Gun  and  Rocket  Propulsion  Systems” ,  Sixteenth  (International) 
Symposium  on  Combustion >  1177-1192  (1976). 

Q 

C.  W.  Nelson "Response  of  Three  Types  of  Transient  Combustion  Models 
to  Gun  Pressurization" 3  Combustion  and  Flame 32j  317-319  (1978). 
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and  Nelson7  showed  that  three  thermal  theory  models  of  the  KTSS8  type 
give  essentially  the  same  response. 

Most  of  these  calculations  assumed  (1)  a  quasi-steady  gas  phase 
with  uniform  heat  release  and  (2)  constant  thermal  properties  of  both 
solid  and  gas.  The  uniform  gas  phase  assumption  derived  from  a  concept 
of  a  diffusion  controlled  flame  of  composite  propellants  not  generally 
held  applicable  to  homogeneous  propellants  with  kinetically  controlled 
gas  phase  reaction.  Thermal  properties  of  the  solid  have  been  found, 
in  the  few  measurements  made,9>10  to  depend  on  temperature.  The  quasi¬ 
steady  gas  phase  remains  the  only  tractable  approach  for  reasonable 
computing  until  methods  like  those  of  Kookerl-  can  be  applied  to  solid 
propellant  combustion  and  considerably  simplified.  Attempts  at  simpli¬ 
fied  unsteady  gas  phase  treatments  tend  to  make  critical  but  indefensible 
assumptions. 12, 13 

14 

Nelson,  et  al.  have  shown  that  predicted  pressure  wave  develop¬ 
ment  in  gun  chambers  is  magnified  by  using  a  Zeldovich  transient  burning 
model  to  replace  the  quasi-steady  regression.  DDT  calculations  (e.g.. 


D.  E.  Kooker  and  C.  W.  Nelson,  "Numerical  Solution  of  Three  Solid 
Propellant  Combustion  Models  During  a  Gun  Pressure  Transient ",  USA 
Ballistic  Research  Laboratory ,  Report  1953  (1977).  (See  also  ASME 
Journal  of  Beat  Transfer,  in  press).  (AD  #A035250) 

Q 

H.  Krier,  J.  S.  T'ien,  W.  A.  Sirignano  and  M.  Summerfield,  "Non-Steady 
Burning  Phenomenon  of  Solid  Propellants:  Theory  and  Experiments",  AIM 
Journal,  6_,  278-285  (1968). 

Q 

J.  R.'  Ward,  "Determination  of  the  Beat  Capacities  of  Gun  Propellants 
ty  Differential  Scanning  Calorimetry ",  Analytical  Calorimetry ,  4, 
143-153  (1977).  “ 

10 

M.  L.  Cohen,  M.  Imber  "Variation  of  the  Thermal  Conductivity  of  Solid 
Propellants  in  an  Acoustic  Field",  AIM  Journal  5,  1295-1299  (1967). 

11 

D.  E.'  Kooker,  "Numerical  Predictions  for  Laminar  Flame  Propagation  in 
Confined  Ozone/Oxygen  Mixtures:  Influence  of  Initial  Temperature  and 
Pressure ",  AIM  Paper  79-0292,  New  Orleans  (1979). 

12 

C.  W.  Nelson,  "Another  Comment  on  the  Transient  Burning  Rate  Model  of 
Suhas  and  Bose",  Combustion  and  Flame,  in  press. 

13 

H,  K.  Suhas  and  T .  K.  Bosey  nA  Mathematical  Model  to  'Predict  Transient 
Burning  Bate  and  Pressure  Decay  Rates  for  Extinction  of  Composite  Pro- 
pellants ",  Combustion  and  Flamey  28 ?  145-153  (1977). 

14 

C.  W.  Nelson 3  P.  S.  Gough 3  and  F .  W.  Robbins >  !t Ignition  Transients  in 
Flame  spreading  Calculations™ y  AIAA  Paper  79—0166 ^  New  Orleans  (1979). 
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Beckstead,  et  al.15)  where  steep  pressure  fronts  are  well  documented 
typically  use  a  quasi-steady  regression  law.  ■ 

The  purpose  of  this  short  report  is  to  examine  theoretically  the 
effect  of  an  assumption  of  variable  thermal  properties  on  the  transient 
regression  rates  of  solid  propellants. 


II.  THEORY 

Thermal  theory  models  solve  the  energy  equation  in  the  solid, 

9T  9T  9  9T-. 

pc  at  +  rpc  ET  "  ■  CU 

7 

which  in  the  nondimensional  variables  of  Kooker  and  Nelson  plus 


L  =  X/X 

o 

C  =  c/c 

o 

becomes 


90  90  _  1  9  ,.90. 

9t  +  R  9^  "  C  9TT  (V  - 


Boundary  conditions  for  the  solid  are: 


C2) 

(3) 

(4) 


0  =  0  as  n  =  -°° 


(5) 


=  f[Ls,  Cs,  R,  PC?)]  at  n  =  0. 

The  specific  form  of  the  surface  boundary  condition  is  set  by  the  choice 
of  gas  phase  model. 

The  initial  condition  is  one  of  steady  state  burning  at  an  initial 
pressure.  In  problems  where  transients  seem  likely  to  be  important, 
the  actual  initial  condition  is  a  cold  propellant.  A  complete  calculation 
would  include  ignition  and  the  transition  to  self-sustaining  combustion* 


15 

M.  W.  Beckstead,  N.  L.  Peterson ,  D.  T.  Pilcher;  and  B.  D.  Hopkins; 
"Convective  Combustion  Modeling  Applied  to  Deflagration  to  Detonation 
Transition  of  HMX";  Combustion  and  Flame  30;  231-241  (1977). 
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Unfortunately,  since  there  is  no  useful  model  of  the  transition,  the 
problem  is  avoided  by  assuming  steady  state  combustion  has  already  been 
achieved. 


At  steady  state,  the  temperature  profile  is  given  by  the  solution 
of  the  energy  equation 


Ri!  -  I  JL  n  39'i 
an  “  c  an  CIW 


(6) 


subject  to  boundary  conditions 


0=0  as  n  4 

0=1  at  n  =  0. 

Note  that  without  the  temperature  dependent  properties,  the  initial 
condition  would  be  the  familiar 

0  =  en 

With  variable  properties,  a  numerical  solution  of  Eq.  (6)  provides  the 
initial  condition. 

In  contrast  with  the  constant  properties  case  where  the  initial 
value  of  the  surface  gradient  is 


the  variable  properties  value  is  approximately 


for  the  specific  values  of  this  problem.  It  will  in  general  vary  from 
unity  in  the  direction  opposing  the  change  in  the  thermal  diffusivity. 

One  simple  treatment  of  the  gas  phase  is  to  assume  a  uniform 
reaction  rate  as  has  been  done  by  several  models  (e.g.,  KTSS8) .  With 
the  reaction  rate  so  specified,  the  gas  phase  energy  equation  can  then 
be  integrated  to  find  the  conduction  term  at  the  solid  boundary.  The 
Kooker-Zinnlb  model  here  yields  the  boundary  condition 


16 

D.  E.  Hooker  and  B.  T.  Zinn,  "Numerical  Investigation  of  Nonlinear 
Axial  Instabilities  in  Solid  Rooket  Motors ",  USA  Ballistic  Research 
Laboratories  Contract  Report  141  (1974).  -(AD  #776954) 
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9±  _  R  rH 

9n  ~  L  tH 


ces-i) 


(C 


c 

E 

c 

0 


)] 


ZP 

RL 


The  constant  Z  is  determined  from  the  initial  condition.  H  is  the  non- 
dimensional  surface  heat  release. 


H  =  Q  /c  (T  -T  ) 
xs  ov  so  0J 

Another  treatment  is  to  ignore  the  gas  phase  details  as  is  done  in 
the  Zeldovich  approach  wherein  it  is  assumed  that  the  functional  form 
of  the  heat ' feedback  is  the  same  for  any  combination  of  regression  rate 
and  pressure. 

For  constant  properties,  the  heat  feedback  term  is  the  transformed 
version  of  the  Zeldovich  boundary  condition, 6 


i£  _  R  rfl  _ L. 

3n  “  L  1  s  "a  (T 

p  so 

For  variable  properties, 
that  at  steady  state 


i  R 
-T  )  ln  R 
oJ  s 


the  term  is 


more  complicated  after 


considering 


C (0)d0  . 


The  essence  of  the  Zeldovich  treatment  is  to  replace  0  with  an  equiva¬ 
lent  0eq  to  give  the  same  regression  rate-pressure  relationship.  For  a 
linear  heat  capacity  dependence  (C  =  Ca  +  CgO) ,  the  surface  gradient 
then  becomes 


P-  =  R[C  (0  -0  )  +  (02-02  ) 

L  Av  s  eqJ  2  v  s  eqJ 


where 


1  .  R 

In  tt— 


eq  a  (T  -T  )  R 
P  s„  oJ 
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III.  MEASURED  THERMAL  PROPERTIES 


Ward  measured  temperature  dependence  of  heat  capacity  below  the 
decomposition  point  for  nitrocellulose  and  X14  propellant.  For  the  X14 
propellant  the  reported  value  was 


C  =  .118  +  .66xl0'3  T (K)  cal/gK 


over  the  temperature  range  283-343  K. 

A  difficulty  enters  here  because  the  theory  has  assumed  an  inert 
solid  at  all  temperatures.  Reactions  in  the  solid  which  prevent  heat 
capacity  measurement  have  been  assumed  to  occur  only  in  a  collapsed  zone 
at  the  surface.  A  simple  solution  is  to  extrapolate  the  measurements 
to  the  calculated  temperature,  however  high.  The  alternative  is  a  com¬ 
plete  recasting  of  the  model  which  would  require  knowledge  of  the 
chemical  reactions  in  the  solid. 

Cohen  measured  the  thermal  conductivity  of  two  double  base  pro¬ 
pellants.  For  DQO  propellant  the  value  is 


A  -  1.90x10  +  2.46x10  ^  T(K)  cal/cmsecK 


over  the  temperature  range  267-317  K.  Again,  the  simplicity  of  an  extra 
polation  outweighs  the  more  difficult  alternatives. 


IV.  RESULTS 

With  the  variable  properties,  transients  are  less  pronounced  than 
with  constant,  low  temperature  properties.  Figure  1  shows  the  effect  of 
property  variation  on  a  calculation  reported  by  Kooker  and  Nelson?  for 
a  distributed  flame.  It  shows  that  transient  response  is  muted  by  the 
variable  properties.  Pressure  was  rising  monotonically  from  7  to  20  MPa 
during  the  transient  as  taken  from  the  records  of  a  105  mm  tank  gun 
which  reaches  400  MPa  in  about  3  ms.  In  Figure  1,  CONCOLD  means  con¬ 
stant  properties  evaluated  at  the  cold  boundary;  CONAVG  means  constant 
properties  evaluated  at  an  interim  temperature  (460  K) ;  VAR  means  the 
variable  properties. 
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KOOKER  MODEL 
Q$  =  82  cal/g 
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2  3  4  5 

PRESSURE  (ND) 


Figure  1.  Effect  of  Variable  Properties 


The  effect  of  pressurization  rate  on  the  Kooker  model  is  shown  in 
Figure  2.  In  the  105  mm  tank  gun,  the  pressurization  rate  at  the  7-20 
MPa  levels  is  about  7000  MPa/sec.  The  pressurization  rate  of  a  grain 
during  the  flamespreading  process  in  a  155  mm  howitzer  is  calculated  to 
be  about  14,000  MPa/sec  at  the  middle  of  the  bed  and  about  40,000  MPa/sec 
at  the  stagnation  of  bed  against  the  projectile  base.  If  the  Kooker 
flame  assumptions  applied,  these  results  suggest  that  burning  rate 
excursions  would  extend  to  pressures  over  100  MPa.  With  the  high 
pressurization  rates  of  DDT,  such  an  effect  cannot  be  ignored  on  the 
assumption  that  the  excursions  apply  only  at  low  pressures. 
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Figure  2.  Effect  of  Pressurization  Rate 


Predictions  of  the  Zeldovich  model  are  shown  in  Figure  3  which 
shows  the  effect  of  varying  temperature  sensitivity  and  pressurization 
rate.  Of  note  is  the  slow  return  of  the  relative  rate  toward  unity. 

For  the  gun  pressurization,  an  excursion  of  over  50%  is  predicted  to 
survive  for  the  whole  ballistic  cycle.  With  a  steep  pressurization 
rate  (7x10-5  MPa/sec)  the  excursion  is  rapid  on  the  time  scale  of  a  gun 
but  does  not  peak  until  50  MPa.  Notable  here  is  the  transition  from  a 
lag  (25%  of  quasi-steady  rate)  to  an  overshoot  of  150%  in  the  entire 
calculated  cycle  to  a  non-dimensional  pressure  of  173  (about  1.2  Kbars) . 
For  DDT  analysis,  these  results  say  that  a  quasi-steady  assumption 
ignores  the  predicted  highly  transient  behavior. 

Surface  heat  release  is  a  critical  parameter  in  the  Kooker  type 
model;  the  non-dimensional  heat  release  is 


H  =  Q  /c 
xs  o 


(T  -T  ) 
s  oJ 
o 


The  effect  of  varying  H  is  shown  in  Figure  4.  The  response  is  similar 
to  that  of  the  constant  properties  solutions  of  Kooker  and  Nelson. 7 
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Experimental  values  for  Qs  are  not  well  known;  Kubota,  et  al.17  have 
suggested  experimental  values  which  vary  with  regression  rate. 


15 


The  model's  response  to  a  step  doubling  of  pressure  from  1  to  2  is 
shown  in  Figure  5.  Such  a  pressure  transient  represents  a  limit  of  a 
steep  pressure  wave.  The  step  used  here  has  no  immediate  practical 
value;  it  is  intended  only  as  a  test  of  the  regression  rate  response 
lag. 

STEP  PRESSURE  RESPONSE 


CONCLUSIONS 

1.  The  regression  rate  overshoots  predicted  with  constant  initial 

temperature  thermal  properties  are  muted  by  the  substitution  of  variable 
thermal  properties. 

2.  Higher  pressurization  rates  lead  to  higher  transients  and  higher 
pressures  at  the  peak  overshoot. 

3.  No  simple  correlation  of  regression  rate  excursion  with  pres¬ 
surization  rate. 

4.  Predicted  excursions  are  dependent  on  assumed  surface  heat 
release. 


5.  Regression  rate  excursions  of  two  times  the  quasi-steady  rate 
are  predicted. 
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C  heat  capacity  (nd)  c/co 

cq  reference  heat  capacity 

0^  gas  phase  heat  capacity 

H  surface  heat  release  (nd) 

CA 

£  constants  in  heat  capacity  temperature  dependence 

B 

L  conductivity  (nd)  X/Xq 

P  pressure  (nd)  P/Pq 

Pq  initial  pressure 

r  regression  rate 

R  regression  rate  (nd)  r/rQ 

r^  steady  state  regression  rate 

rQ  initial  regression  rate 

T  temperature 

Tq  reference  temperature 

reference  surface  temperature 

aQ  reference  thermal  diffusivity 

n  distance  (nd)  xr  /a 

X  thermal  conductivity 

reference  thermal  conductivity 

p  density 

0  temperature  (nd)  CT-Tq) /tS0“t03 

Op  temperature  sensitivity  of  regression  rate 

2 

t  time  (nd)  t  rQ/a 
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